PERCEPTION predicts patient response and resistance to treatment using single-cell transcriptomics of their tumors.

Tailoring optimal treatment for individual cancer patients remains a significant challenge. To address this issue, we developed PERCEPTION (PERsonalized Single-Cell Expression-Based Planning for Treatments In ONcology), a precision oncology computational pipeline. Our approach uses publicly available matched bulk and single-cell (sc) expression profiles from large-scale cell-line drug screens. These profiles help build treatment response models based on patients' sc-tumor transcriptomics. PERCEPTION demonstrates success in predicting responses to targeted therapies in cultured and patient-tumor-derived primary cells, as well as in two clinical trials for multiple myeloma and breast cancer. It also captures the resistance development in patients with lung cancer treated with tyrosine kinase inhibitors. PERCEPTION outperforms published state-of-the-art sc-based and bulk-based predictors in all clinical cohorts. PERCEPTION is accessible at https://github.com/ruppinlab/PERCEPTION . Our work, showcasing patient stratification using sc-expression profiles of their tumors, will encourage the adoption of sc-omics profiling in clinical settings, enhancing precision oncology tools based on sc-omics.

Charting the transcriptomic landscape of primary and metastatic cancers in relation to their origin and target normal tissues.

Metastasis is a leading cause of cancer-related deaths, yet understanding how metastatic tumors adapt from their origin to target tissues is challenging. To address this, we assessed whether primary and metastatic tumors resemble their tissue of origin or target tissue in terms of gene expression. We analyzed gene expression profiles from various cancer types, including single-cell and bulk RNA-seq data, in both paired and unpaired primary and metastatic patient cohorts. We quantified the transcriptomic distances between tumor samples and their normal tissues, revealing that primary tumors are more similar to their tissue of origin, while metastases shift towards the target tissue. Pathway-level analysis highlighted critical transcriptomic changes during metastasis. Notably, primary cancers exhibited higher activity in cancer hallmarks, including Activating Invasion and Metastasis , compared to metastatic cancers. This comprehensive landscape analysis provides insight into how cancer tumors adapt to their metastatic environments, providing a transcriptome-wide view of the processes involved.

Prediction of cancer treatment response from histopathology images through imputed transcriptomics.

Advances in artificial intelligence have paved the way for leveraging hematoxylin and eosin (H&E)-stained tumor slides for precision oncology. We present ENLIGHT-DeepPT, an approach for predicting response to multiple targeted and immunotherapies from H&E-slides. In difference from existing approaches that aim to predict treatment response directly from the slides, ENLIGHT-DeepPT is an indirect two-step approach consisting of (1) DeepPT, a new deep-learning framework that predicts genome-wide tumor mRNA expression from slides, and (2) ENLIGHT, which predicts response based on the DeepPT inferred expression values. DeepPT successfully predicts transcriptomics in all 16 TCGA cohorts tested and generalizes well to two independent datasets. Our key contribution is showing that ENLIGHT-DeepPT successfully predicts true responders in five independent patients' cohorts involving four different treatments spanning six cancer types with an overall odds ratio of 2.44, increasing the baseline response rate by 43.47% among predicted responders, without the need for any treatment data for training. Furthermore, its prediction accuracy on these datasets is comparable to a supervised approach predicting the response directly from the images, which needs to be trained and tested on the same cohort. ENLIGHT-DeepPT future application could provide clinicians with rapid treatment recommendations to an array of different therapies and importantly, may contribute to advancing precision oncology in developing countries.

Pan-Cancer Analysis of Patient Tumor Single-Cell Transcriptomes Identifies Promising Selective and Safe Chimeric Antigen Receptor Targets in Head and Neck Cancer.

Chimeric antigen receptor (CAR) T cell therapies have yielded transformative clinical successes for patients with blood tumors, but their full potential remains to be unleashed against solid tumors. One challenge is finding selective targets, which we define intuitively to be cell surface proteins that are expressed widely by cancer cells but minimally by healthy cells in the tumor microenvironment and other normal tissues. Analyzing patient tumor single-cell transcriptomics data, we first defined and quantified selectivity and safety scores of existing CAR targets for indications in which they are in clinical trials or approved. We then sought new candidate cell surface CAR targets that have better selectivity and safety scores than those currently being tested. Remarkably, in almost all cancer types, we could not find such better targets, testifying to the near optimality of the current target space. However, in human papillomavirus (HPV)-negative head and neck squamous cell carcinoma (HNSC), for which there is currently a dearth of existing CAR targets, we identified a total of twenty candidate novel CAR targets, five of which have both superior selectivity and safety scores. These newly identified cell surface targets lay a basis for future investigations that may lead to better CAR treatments in HNSC.

The GPCR-Gαs-PKA signaling axis promotes T cell dysfunction and cancer immunotherapy failure.

Immune checkpoint blockade (ICB) targeting PD-1 and CTLA-4 has revolutionized cancer treatment. However, many cancers do not respond to ICB, prompting the search for additional strategies to achieve durable responses. G-protein-coupled receptors (GPCRs) are the most intensively studied drug targets but are underexplored in immuno-oncology. Here, we cross-integrated large singe-cell RNA-sequencing datasets from CD8+ T cells covering 19 distinct cancer types and identified an enrichment of Gαs-coupled GPCRs on exhausted CD8+ T cells. These include EP2, EP4, A2AR, ß1AR and ß2AR, all of which promote T cell dysfunction. We also developed transgenic mice expressing a chemogenetic CD8-restricted Gαs-DREADD to activate CD8-restricted Gαs signaling and show that a Gαs-PKA signaling axis promotes CD8+ T cell dysfunction and immunotherapy failure. These data indicate that Gαs-GPCRs are druggable immune checkpoints that might be targeted to enhance the response to ICB immunotherapies.

Predicting tumor immune microenvironment and checkpoint therapy response of head & neck cancer patients from blood immune single-cell transcriptomics.

The immune state of tumor microenvironment is crucial for determining immunotherapy response but is not readily accessible. Here we investigate if we can infer the tumor immune state from the blood and further predict immunotherapy response. First, we analyze a dataset of head and neck squamous cell carcinoma (HNSCC) patients with matched scRNA-Seq of peripheral blood mononuclear cells (PBMCs) and tumor tissues. We find that the tumor immune cell fractions of different immune cell types and many of the genes they express can be inferred from the matched PBMC scRNA-Seq. Second, analyzing another HNSCC dataset with PBMC scRNA-Seq and immunotherapy response, we find that the inferred ratio between tumor memory B and regulatory T cell fractions is predictive of immunotherapy response and is superior to the well-established cytolytic and exhausted T-cell signatures. Overall, these results showcase the potential of scRNA-Seq liquid biopsies in cancer immunotherapy, calling for their larger-scale testing. Significance: This head and neck cancer study demonstrates the potential of using blood single-cell transcriptomics to (1) infer the tumor immune status and (2) predict immunotherapy response from the tumor immune status inferred from blood. These results showcase the potential of single-cell transcriptomics liquid biopsies for further advancing personalized cancer immunotherapy.

Clinically oriented prediction of patient response to targeted and immunotherapies from the tumor transcriptome.

BACKGROUND: Precision oncology is gradually advancing into mainstream clinical practice, demonstrating significant survival benefits. However, eligibility and response rates remain limited in many cases, calling for better predictive biomarkers. METHODS: We present ENLIGHT, a transcriptomics-based computational approach that identifies clinically relevant genetic interactions and uses them to predict a patient's response to a variety of therapies in multiple cancer types without training on previous treatment response data. We study ENLIGHT in two translationally oriented scenarios: personalized oncology (PO), aimed at prioritizing treatments for a single patient, and clinical trial design (CTD), selecting the most likely responders in a patient cohort. FINDINGS: Evaluating ENLIGHT's performance on 21 blinded clinical trial datasets in the PO setting, we show that it can effectively predict a patient's treatment response across multiple therapies and cancer types. Its prediction accuracy is better than previously published transcriptomics-based signatures and is comparable with that of supervised predictors developed for specific indications and drugs. In combination with the interferon-γ signature, ENLIGHT achieves an odds ratio larger than 4 in predicting response to immune checkpoint therapy. In the CTD scenario, ENLIGHT can potentially enhance clinical trial success for immunotherapies and other monoclonal antibodies by excluding non-responders while overall achieving more than 90% of the response rate attainable under an optimal exclusion strategy. CONCLUSIONS: ENLIGHT demonstrably enhances the ability to predict therapeutic response across multiple cancer types from the bulk tumor transcriptome. FUNDING: This research was supported in part by the Intramural Research Program, NIH and by the Israeli Innovation Authority.

Sex Biases in Cancer and Autoimmune Disease Incidence Are Strongly Positively Correlated with Mitochondrial Gene Expression across Human Tissues.

Cancer occurs more frequently in men while autoimmune diseases (AIDs) occur more frequently in women. To explore whether these sex biases have a common basis, we collected 167 AID incidence studies from many countries for tissues that have both a cancer type and an AID that arise from that tissue. Analyzing a total of 182 country-specific, tissue-matched cancer-AID incidence rate sex bias data pairs, we find that, indeed, the sex biases observed in the incidence of AIDs and cancers that occur in the same tissue are positively correlated across human tissues. The common key factor whose levels across human tissues are most strongly associated with these incidence rate sex biases is the sex bias in the expression of the 37 genes encoded in the mitochondrial genome.

Dysregulation of HNF1B/Clusterin axis enhances disease progression in a highly aggressive subset of pancreatic cancer patients.

Pancreatic ductal adenocarcinoma (PDAC) is a lethal malignancy and is largely refractory to available treatments. Identifying key pathways associated with disease aggressiveness and therapeutic resistance may characterize candidate targets to improve patient outcomes. We used a strategy of examining the tumors from a subset of PDAC patient cohorts with the worst survival to understand the underlying mechanisms of aggressive disease progression and to identify candidate molecular targets with potential therapeutic significance. Non-negative matrix factorization (NMF) clustering, using gene expression profile, revealed three patient subsets. A 142-gene signature specific to the subset with the worst patient survival, predicted prognosis and stratified patients with significantly different survival in the test and validation cohorts. Gene-network and pathway analysis of the 142-gene signature revealed dysregulation of Clusterin (CLU) in the most aggressive patient subset in our patient cohort. Hepatocyte nuclear factor 1 b (HNF1B) positively regulated CLU, and a lower expression of HNF1B and CLU was associated with poor patient survival. Mechanistic and functional analyses revealed that CLU inhibits proliferation, 3D spheroid growth, invasiveness and epithelial-to-mesenchymal transition (EMT) in pancreatic cancer cell lines. Mechanistically, CLU enhanced proteasomal degradation of EMT-regulator, ZEB1. In addition, orthotopic transplant of CLU-expressing pancreatic cancer cells reduced tumor growth in mice. Furthermore, CLU enhanced sensitivity of pancreatic cancer cells representing aggressive patient subset, to the chemotherapeutic drug gemcitabine. Taken together, HNF1B/CLU axis negatively regulates pancreatic cancer progression and may potentially be useful in designing novel strategies to attenuate disease progression in PDAC patients.

Big data in basic and translational cancer research.

Historically, the primary focus of cancer research has been molecular and clinical studies of a few essential pathways and genes. Recent years have seen the rapid accumulation of large-scale cancer omics data catalysed by breakthroughs in high-throughput technologies. This fast data growth has given rise to an evolving concept of 'big data' in cancer, whose analysis demands large computational resources and can potentially bring novel insights into essential questions. Indeed, the combination of big data, bioinformatics and artificial intelligence has led to notable advances in our basic understanding of cancer biology and to translational advancements. Further advances will require a concerted effort among data scientists, clinicians, biologists and policymakers. Here, we review the current state of the art and future challenges for harnessing big data to advance cancer research and treatment.

MCP-1/CCR2 axis inhibition sensitizes the brain microenvironment against melanoma brain metastasis progression.

Development of resistance to chemo- and immunotherapies often occurs following treatment of melanoma brain metastasis (MBM). The brain microenvironment (BME), particularly astrocytes, cooperate toward MBM progression by upregulating secreted factors, among which we found that monocyte chemoattractant protein-1 (MCP-1) and its receptors, CCR2 and CCR4, were overexpressed in MBM compared with primary lesions. Among other sources of MCP-1 in the brain, we show that melanoma cells altered astrocyte secretome and evoked MCP-1 expression and secretion, which in turn induced CCR2 expression in melanoma cells, enhancing in vitro tumorigenic properties, such as proliferation, migration, and invasion of melanoma cells. In vivo pharmacological blockade of MCP-1 or molecular knockout of CCR2/CCR4 increased the infiltration of cytotoxic CD8+ T cells and attenuated the immunosuppressive phenotype of the BME as shown by decreased infiltration of Tregs and tumor-associated macrophages/microglia in several models of intracranially injected MBM. These in vivo strategies led to decreased MBM outgrowth and prolonged the overall survival of the mice. Our findings highlight the therapeutic potential of inhibiting interactions between BME and melanoma cells for the treatment of this disease.

AXL and Error-Prone DNA Replication Confer Drug Resistance and Offer Strategies to Treat EGFR-Mutant Lung Cancer.

Anticancer therapies have been limited by the emergence of mutations and other adaptations. In bacteria, antibiotics activate the SOS response, which mobilizes error-prone factors that allow for continuous replication at the cost of mutagenesis. We investigated whether the treatment of lung cancer with EGFR inhibitors (EGFRi) similarly engages hypermutators. In cycling drug-tolerant persister (DTP) cells and in EGFRi-treated patients presenting residual disease, we observed upregulation of GAS6, whereas ablation of GAS6's receptor, AXL, eradicated resistance. Reciprocally, AXL overexpression enhanced DTP survival and accelerated the emergence of T790M, an EGFR mutation typical to resistant cells. Mechanistically, AXL induces low-fidelity DNA polymerases and activates their organizer, RAD18, by promoting neddylation. Metabolomics uncovered another hypermutator, AXL-driven activation of MYC, and increased purine synthesis that is unbalanced by pyrimidines. Aligning anti-AXL combination treatments with the transition from DTPs to resistant cells cured patient-derived xenografts. Hence, similar to bacteria, tumors tolerate therapy by engaging pharmacologically targetable endogenous mutators. SIGNIFICANCE: EGFR-mutant lung cancers treated with kinase inhibitors often evolve resistance due to secondary mutations. We report that in similarity to the bacterial SOS response stimulated by antibiotics, endogenous mutators are activated in drug-treated cells, and this heralds tolerance. Blocking the process prevented resistance in xenograft models, which offers new treatment strategies. This article is highlighted in the In This Issue feature, p. 2483.

Membrane-Tethered Mucin 1 Is Stimulated by Interferon and Virus Infection in Multiple Cell Types and Inhibits Influenza A Virus Infection in Human Airway Epithelium.

Influenza A virus (IAV) causes significant morbidity and mortality in the human population. Tethered mucin 1 (MUC1) is highly expressed in airway epithelium, the primary site of IAV replication, and also by other cell types that influence IAV infection, including macrophages. MUC1 has the potential to influence infection dynamics through physical interactions and/or signaling activity, yet MUC1 modulation and its impact during viral pathogenesis remain unclear. Thus, we investigated MUC1-IAV interactions in an in vitro model of human airway epithelium (HAE). Our data indicate that a recombinant IAV hemagglutinin (H3) and H3N2 virus can bind endogenous HAE MUC1. Notably, infection of HAE with H1N1 or H3N2 IAV strains does not trigger MUC1 shedding but instead stimulates an increase in cell-associated MUC1 protein. We observed a similar increase after type I or III interferon (IFN) stimulation; however, inhibition of IFN signaling during H1N1 infection only partially abrogated this increase, indicating that multiple soluble factors contribute to MUC1 upregulation during the antiviral response. In addition to HAE, primary human monocyte-derived macrophages also upregulated MUC1 protein in response to IFN treatment and conditioned media from IAV-infected HAE. Then, to determine the impact of MUC1 on IAV pathogenesis, we developed HAE genetically depleted of MUC1 and found that MUC1 knockout cultures exhibited enhanced viral growth compared to control cultures for several IAV strains. Together, our data support a model whereby MUC1 inhibits productive uptake of IAV in HAE. Infection then stimulates MUC1 expression on multiple cell types through IFN-dependent and -independent mechanisms that further impact infection dynamics. IMPORTANCE Influenza A virus (IAV) targets airway epithelial cells for infection. Large, heavily glycosylated molecules known as tethered mucins extend from the airway epithelial cell surface and may physically restrict pathogen access to underlying cells. Additionally, tethered mucin 1 (MUC1) can be differentially phosphorylated based on external stimuli and can influence inflammation. Given MUC1's multifunctional capability, we sought to define its role during IAV infection. Here, we demonstrate that IAV directly interacts with MUC1 in a physiologically relevant model of human airway epithelium (HAE) and find that MUC1 protein expression is elevated throughout the epithelium and in primary human monocyte-derived macrophages in response to antiviral signals produced during infection. Using CRISPR/Cas9-modified HAE, we demonstrated more efficient IAV infection when MUC1 is genetically ablated. Our data suggest that MUC1 physically restricts IAV uptake and represents a dynamic component of the host response that acts to inhibit viral spread, yielding new insight into mucin-mediated antiviral defense.

Synthetic lethality-based prediction of anti-SARS-CoV-2 targets.

Novel strategies are needed to identify drug targets and treatments for the COVID-19 pandemic. The altered gene expression of virus-infected host cells provides an opportunity to specifically inhibit viral propagation via targeting the synthetic lethal and synthetic dosage lethal (SL/SDL) partners of such altered host genes. Pursuing this disparate antiviral strategy, here we comprehensively analyzed multiple in vitro and in vivo bulk and single-cell RNA-sequencing datasets of SARS-CoV-2 infection to predict clinically relevant candidate antiviral targets that are SL/SDL with altered host genes. The predicted SL/SDL-based targets are highly enriched for infected cell inhibiting genes reported in four SARS-CoV-2 CRISPR-Cas9 genome-wide genetic screens. We further selected a focused subset of 26 genes that we experimentally tested in a targeted siRNA screen using human Caco-2 cells. Notably, as predicted, knocking down these targets reduced viral replication and cell viability only under the infected condition without harming noninfected healthy cells.

Immune Determinants of the Association between Tumor Mutational Burden and Immunotherapy Response across Cancer Types.

The FDA has recently approved a high tumor mutational burden (TMB-high) biomarker, defined by ≥10 mutations/Mb, for the treatment of solid tumors with pembrolizumab, an immune checkpoint inhibitor (ICI) that targets PD1. However, recent studies have shown that this TMB-high biomarker is only able to stratify ICI responders in a subset of cancer types, and the mechanisms underlying this observation have remained unknown. The tumor immune microenvironment (TME) may modulate the stratification power of TMB (termed TMB power), determining if it will be predictive of ICI response in a given cancer type. To systematically study this hypothesis, we inferred the levels of 31 immune-related factors characteristic of the TME of different cancer types in The Cancer Genome Atlas. Integration of this information with TMB and response data of 2,277 patients treated with anti-PD1 identified key immune factors that determine TMB power across 14 different cancer types. We find that high levels of M1 macrophages and low resting dendritic cells in the TME characterized cancer types with high TMB power. A model based on these two immune factors strongly predicted TMB power in a given cancer type during cross-validation and testing (Spearman Rho = 0.76 and 1, respectively). Using this model, we predicted the TMB power in nine additional cancer types, including rare cancers, for which TMB and ICI response data are not yet publicly available. Our analysis indicates that TMB-high may be highly predictive of ICI response in cervical squamous cell carcinoma, suggesting that such a study should be prioritized. SIGNIFICANCE: This study uncovers immune-related factors that may modulate the relationship between high tumor mutational burden and ICI response, which can help prioritize cancer types for clinical trials.

The landscape of receptor-mediated precision cancer combination therapy via a single-cell perspective.

Mining a large cohort of single-cell transcriptomics data, here we employ combinatorial optimization techniques to chart the landscape of optimal combination therapies in cancer. We assume that each individual therapy can target any one of 1269 genes encoding cell surface receptors, which may be targets of CAR-T, conjugated antibodies or coated nanoparticle therapies. We find that in most cancer types, personalized combinations composed of at most four targets are then sufficient for killing at least 80% of tumor cells while sparing at least 90% of nontumor cells in the tumor microenvironment. However, as more stringent and selective killing is required, the number of targets needed rises rapidly. Emerging individual targets include PTPRZ1 for brain and head and neck cancers and EGFR in multiple tumor types. In sum, this study provides a computational estimate of the identity and number of targets needed in combination to target cancers selectively and precisely.

A systematic genome-wide mapping of oncogenic mutation selection during CRISPR-Cas9 genome editing.

Recent studies have reported that genome editing by CRISPR-Cas9 induces a DNA damage response mediated by p53 in primary cells hampering their growth. This could lead to a selection of cells with pre-existing p53 mutations. In this study, employing an integrated computational and experimental framework, we systematically investigated the possibility of selection of additional cancer driver mutations during CRISPR-Cas9 gene editing. We first confirm the previous findings of the selection for pre-existing p53 mutations by CRISPR-Cas9. We next demonstrate that similar to p53, wildtype KRAS may also hamper the growth of Cas9-edited cells, potentially conferring a selective advantage to pre-existing KRAS-mutant cells. These selective effects are widespread, extending across cell-types and methods of CRISPR-Cas9 delivery and the strength of selection depends on the sgRNA sequence and the gene being edited. The selection for pre-existing p53 or KRAS mutations may confound CRISPR-Cas9 screens in cancer cells and more importantly, calls for monitoring patients undergoing CRISPR-Cas9-based editing for clinical therapeutics for pre-existing p53 and KRAS mutations.

Genome-scale metabolic modeling reveals SARS-CoV-2-induced metabolic changes and antiviral targets.

Tremendous progress has been made to control the COVID-19 pandemic caused by the SARS-CoV-2 virus. However, effective therapeutic options are still rare. Drug repurposing and combination represent practical strategies to address this urgent unmet medical need. Viruses, including coronaviruses, are known to hijack host metabolism to facilitate viral proliferation, making targeting host metabolism a promising antiviral approach. Here, we describe an integrated analysis of 12 published in vitro and human patient gene expression datasets on SARS-CoV-2 infection using genome-scale metabolic modeling (GEM), revealing complicated host metabolism reprogramming during SARS-CoV-2 infection. We next applied the GEM-based metabolic transformation algorithm to predict anti-SARS-CoV-2 targets that counteract the virus-induced metabolic changes. We successfully validated these targets using published drug and genetic screen data and by performing an siRNA assay in Caco-2 cells. Further generating and analyzing RNA-sequencing data of remdesivir-treated Vero E6 cell samples, we predicted metabolic targets acting in combination with remdesivir, an approved anti-SARS-CoV-2 drug. Our study provides clinical data-supported candidate anti-SARS-CoV-2 targets for future evaluation, demonstrating host metabolism targeting as a promising antiviral strategy.